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ABSTRACT 

In this paper, the estimation of R=P(Y < X), when X and Y are two generalized inverted exponential random 
variables with different parameters is considered. This problem arises naturally in the area of reliability for a 
system with strength X and stress Y. The estimation is made using simple random sampling (SRS) and ranked set 
sampling (RSS) approaches. The maximum likelihood estimator (MLE) of R is derived using both approaches. 
Assuming that the common scale parameter is known, MLEs of R are obtained. Monte Carlo simulations are 
performed to compare the estimators obtained using both approaches. . The properties of these estimators are 
investigated and compared with known estimators based on simple random sample (SRS) data. The comparison 
is based on biases, mean squared errors (MSEs) and the efficiency of the estimators of R based on RSS with 
respect to those based on SRS. The estimators based on RSS is found to dominate those based on SRS. 

KEYWORDS: generalized exponential distribution; reliability; stress-strength; ranked set sampling, simple 
random sampling; maximum likelihood estimators. 

I. Introduction 

The estimation of reliability is a very common problem in statistical literature. The most widely 
approach applied for reliability estimation is the well-known stress-strength model. This model is 
used in many applications of physics and engineering such as strength failure and the system collapse. 
In the stress-strength modeling, R-P(Y < X) is a measure of component reliability when it is 
subjected to random stress Y and has strength X. In this context, R can be considered as a measure of 
system performance and it is naturally arise in electrical and electronic systems. Another 
interpretation can be that, the reliability of the system is the probability that the system is strong 
enough to overcome the stress imposed on it. It may be mentioned that R is of greater interest than 
just reliability since it provides a general measure of the difference between two populations and has 
applications in many areas. For example, if X is the response for a control group, and Y refers to a 
treatment group, R is a measure of the effect of the treatment. In addition, it may be mentioned that R 
equals the area under the receiver operating characteristic (ROC) curve for diagnostic test or 
biomarkers with continuous outcome, see; Bamber, [1]. The ROC curve is widely used, in biological, 
medical and health service research, to evaluate the ability of diagnostic tests or biomarkers and to 
distinguish between two groups of subjects, usually non-diseased and diseased subjects. For complete 
review and more applications of R; see [2-14]. 

Ranked set sampling (RSS) is a sampling protocol that can often be used to improve the cost and 
efficiency for experiments [15]. It is often used when a ranking of the sampling units can be obtained 
cheaply without having to actually measure the characteristics of interest, which may be time 
consuming or costly [16,17]. Such a technique is well received and widely applicable in 
environmental applications, reliability and quality control experiments [18-20]. A modification of 
ranked set sampling (RSS) called moving extremes ranked set sampling (MERSS) was considered for 
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the estimation of the scale parameter of scale distributions [21] and an improved RSS estimator for 
the population mean was obtained [22]. On the other hand, Ozturk has developed two sampling 
designs to create artificially stratified samples using RSS [23] 

Recently, many authors have been interested in estimating R using RSS. For example, Sengupta and 
Mukhuti [24], considered an unbiased estimation of R using RSS for exponential populations. Muttlak 
and co-authors [25], proposed three estimators of R using RSS when X and Y independent one- 
parameter exponential populations. In a RSS procedure, m independent sets of SRS each of size m are 
drawn from the distribution under consideration, these samples are ranked by some auxiliary criterion 
that does not require actual measurements and only the /th smallest observation is quantified from the 
/th set, i = l,2,...,m. This completes a cycle of the sampling. Then, the cycle is repeated k times to 
obtain a ranked set sample of size n-mk . 

In this article, we consider the estimation of R, when X and Y are independent but not identically 
distributed generalized inverted exponential distribution (GIED) variables. Abouammoh and 
Alshingiti [26], introduced a shape parameter in the inverted exponential distribution (IED) to obtain 
the generalized inverted exponential distribution (GIED). They derived many distributional properties 
and reliability characteristics of GIED. Assuming it to be a good lifetime model, they obtained 
maximum likelihood estimators, least square estimators and confidence intervals of the two 
parameters involved. The GIED has the following cumulative distribution function (cdf) and 
probability density function (pdf) for X > 0: 

F(x) = l-[l-exp(-A/x)] a , x>0, (1) 

with 

f(x) = a— exp(-X/x)[l-exp(-X/x)] a ~\ x>0 , (2) 

x 2 

where X > is the scale parameter and a > is the shape parameter. 

The rest of paper is organized as follows: the derivation of the stress-strength (R = P(Y<X)) model for 
the generalized inverted exponential distribution is discussed. In sections 3 and 4, the MLEs of R 
using SRS and RSS approaches are derived.. In section 5 simulation studies and comparison between 
estimators obtained using both approaches are discussed. Finally, the paper is concluded in section 6. 

II. The Stress Strength Model 

Let X and Y be two independent GIED random variables with parameters (X,a) and (X,/3) 

respectively. Therefore, the reliability of the system is given by 
l 

R = P(Y <X) = \f x (x)P(Y <X\X = x)dx 
o 

= ja^- exp(-^/x)[l-exp(-/l/x)] a " 1 (l-[l-exp(-/l/x)]^)Jx, 
x 

= 1 _^ = ^_. (3) 
a + jB a + /3 

III. Maximum Likelihood Estimation of R Using SRS 

Let (X l ,X 2 ,...,X n ) and (Y l ,Y 2 ,...,Y m ) be two independent random samples from GIED (X, a) and 
GIED (X,f3) respectively. The likelihood function of X, a and /3 for the observed samples is 

n n n 

L s (data;X,a,P)=X n a n [T\xf] 1 exp(- E II x() T\[l -exp(-A /x,)] a 
i=l i'=l (=1 

m o i m m ft 1 

x rj} m \_\\ yjr 1 exp(- Z A/yj) rni-expC-l/^)]^- 1 (4) 
7=1 J=l j=l 

Therefore, the log-likelihood function of X, a and /? will be 
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n n j n 

\ogL s =(/« + «)log/l+«loga + mlog/? -2 ^logx,- - £ — + (a-l) J^\og[l-exp(-A/ xi)] 

i=\ i=l x i i=\ 

m m 2 m 

-2£log^-£— +Off-l)£log[l-exp(-A/y J )]. (5) 

The estimators 2 ? , a s and S of the parameters X, a and /3 respectively can be then obtained as 
the solution of the likelihood equations 

<^-±±-±±Ha-l)± eXP( ' 1/X ' ) HP-D± eXP( ' A,yj) = (6) 
X i=l Xi j=l yj j =1 x t [l-exp(-/L/x ( )] j=l yj [l-expC-A/^)] 

-+Ylog[l-exp(-A/x ( )] = (7) 

« i=i 

m 

-+£log[l-exp(-A/y,)] = (8) 

" 7=1 

Once the nonlinear equations (6-8) have been solved with respect to A, a and /3 , the MLEs d s and 
P s will be given as 

n 



(9) 



Xlogtl-expC-^/x,)] 

f=i 

and 

p, = - 5 . do) 

Ilogll-expM,/^)] 

7=1 

where I s is the solution of the nonlinear equation, 

(n + m) £ 1 £ 1 ( w +1)5 exp(-4/x i -) 



i=l 

- ( - + i)f; ct p^^- =0 . (id 

f>g[l-ex P (-Vy 7 .)] My;[l-ex P (-^/y ; )] 

.7=1 

By the invariance property of the MLEs, the MLE of R becomes 

[m JJ log(l - exp(-/L s / x t )] 
= . (12) 

s Ps [m X log(l - exp(-^ /**)+!>! log[l - exp(-l s / y ■ )] 

1=1 7=1 

Assuming the scale parameter is known, the MLEs a s , x , P Sy x and R s x of the parameters a , j5 
and are respectively 

a sX = ^ , (13) 

Xlog[l-exp(-/l/x,-)] 
i=l 



^ , (14) 



m 

£log[l-expM/> y -)] 

7=1 

and 



n 

[mJJ log(l-exp(-/l/x,)] 

(15) 



«*,A=— 

[m log(l - exp(-2 / x t ) + n log[l - exp(-A / yj )] 
/=1 ;=1 
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IV. Maximum Likelihood Estimation Of R Using Rss 

Assume that X^, m i = l,...,m x , j = l,...,r x is a ranked set samples with sample size n x -m x r x 

drawn from GIED with shape and scale parameters a and X respectively. Let 
Y (k m )i >k = l,...,m y ,l =l,...,r y be an independent ranked set sample with sample size n y -m y r v 

drawn from GIED with shape and scale parameters P and X respectively. m x and r x are the set size 
and the number of cycles of the RSS sample drawn from GIED (k, a) , while m and r are the set 
size and the number of cycles of the RSS drawn from GIED (A,/?). For simplification purposes, 
X(t m )j and Y^ k:m yi will be denoted as X y and Y kl respectively. The pdf of the random variables 
and Y kl are given by 

*i(*#> = -a^p(-Z/ Xij )[l - «p(-A/*# ]*<»> (1 - [1 - exp(-A/*, f )"*-' , (16) 

-!)!('«* -0 x| 

(*-l)!(m y -k) y L kl 

where k,cc,/3>0, Xg > , Y kl > 0, j = 1,...,£ and / = L..../J, . The likelihood function of k, a and f3 
for the observed samples is given by 

r x m x 

L r (data;a,k) = K^[ riexp(-/l/Xy)[l-exp(-/l/Xy] a(rajc "' +1) " 1 (l-[l-exp(-/l/Xy] 0: ) m ^"' 
7=1 '=1 

* ft ii P \ exp(-^ / y kl ) [1 - exp(-A / >•*, )/ (m > ~* (1 - [1 - exp(-^ / y H )/) m > "* ( 1 8 ) 
l=i k=i y kI 

Therefore, the log -likelihood function of k, a and (3 will be 

* r x m x r y m y 

LogL r =K + m x r x {\oga+\ogk)+m y r y (log/? + log A)- £ Z^/*//-Z I kly^i 

j=\ i=\ i=i k=l 

r x m \ r y m y 

+ (a(m x -i + 1) -1) Z Zlog[l - «p(-A/^)] + (£(m„ -k + 1) -1)2 Zlog[l - enp(-A/y w )] 
/=i i=i ;=u=i 

+ k - o£ Z 1 °g( 1 - p - ex p(- A / ^ )]" )+ («, - *)£ Z ^ - p - ex p(-^ lykdf) ( 1 9) 

j=l i=l Z=l i=l 

where A" is a constant. This implies that 

m x r x , .'A^ n , r m* [1 - exp(-^ / )f log[l - exp(-^ / Xjj )] (20) 

-=^ + (m r -! + l)2 21og[l-exp(-A,/x,,)] + (m v -0Z Z " s " — = ° V ' 

a ;=l/=l ;=l«=l l-[l-exp(-A. /*,-,■)] 

ra v r >' , , ^n^, n .1/ \-i , / ,.^^[l-exp(-X/.v H )]Plog[l-exp(-X/.y H )] ( 2 1) 

— — + (m Y -k + l)2 Ilog[l-exp(-X./y H )] + (m v -/:)Z Z 5 — = 0' 

P »=i*=i i=i k=\ l-[l-exp(-A./ y kl )f 

m x r x +m v r v r x m x j r y m y r m e xp(-A/x,-,) 

} } -I Z— -I I— + (a(m,-i+l)-l)Z L- 7 



and 



j=li=l x ij l=\k=\yu y =i/=Hy[l-exp(-/l/x (7 ] 

expt-A/y^) ' ^expC-A/xyXl-expC-zl/xy)]"" 1 



+ (y9(m y -it + l)-l)i £ ^ + a(m x -i)Y, Z 

/=l*r=iy«[l-exp(-^/yjfe/)] y=l i=l JC !7 (l-[l-exp(-/l/x, 7 )] a ) 

+ PK-^Z exp(-X/y t/ )[l-exp(->./ yt< )f- 1 = Q . (22 ) 
5 j=l*=l y w (l-[l-expH./y iH )] 15 ) 

To derive the MLEs k r Sc r and f3 r of the parameters k, a and /5 , we have to solve the nonlinear 
equations (20-22) with respect to k, a and [3 . Therefore, the MLE of R is given by 

R r = - ^ r . (23) 

a r + fly 
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Once again, assuming that the scale parameter A is known, the MLEs d rX and fi r ^ of the 

parameters a and /3 are the solution of the likelihood equations 

m x r x . .^a'U. _ , . , v ,. , ^^[l-exp(-X/^.)] drA log[l-exp(-X/^)] (24) 

-^ L ^ + (m r -; + l)I 21og[l-exp(-A,/X;:)] + (m x -0Z Z ^ ^ = v > 

«a ,/=i i=i y=i «=i l-[l-exp(-A,/jt g )] ' 

3-5- + (m, + W Ilog[l-exp(-X/ yH )] + K " ^-^^ lo g [l-exp(-X/ yH)] =Q (25) 

P,-,x i=i*=i l-[l-exp(-X,/y H )f' A 

Thus, the MLE of 7? is given by 

Kx = : ^ (26) 

a r,A + Pr,l 

It is clear that the MLEs of X and a denoted by A r and a r cannot be obtained in closed form. Thus, 

we will use simulation to study the properties of these estimators and the properties of the estimators 
of/?. 

V. Simulation Study 

This section presents a numerical comparison between the MLEs of /?=P(Y < X) using both SRS and 
RSS approaches. The estimation is made through biases and mean squared errors (MSE) of the 

estimators R s x and R r x and the efficiency of the estimator R r x with respect the estimator R s x 
where the efficiency of a parameter 9 2 with respect to the parameter 9 l is given by 

tf(U) = «^ (27) 
MSE(d 2 ) 

The simulation studies are made using MATHEMATICA v. 8 for several combinations of the 
parameters n, m , m x , r x , m , r y and R and for r x = r y =5. The random samples of GIED(a,^) 

and GIED(P,^) are generated using the forms 

x= — - — , jt>0, A,a>0 (28) 

log[l- M 1/a ] 

and 

y= ~~mi~' x ^°> ^A >0 (29) 

logtl-v 1 ^] 

where < u < 1 and < v < 1 are uniform random variables. Different simulations are based on 1000 
replications. The results are shown in Tables 1-3. 

Based on Tables 1-3, one can see that the biases and MSEs of both R s A and R r x decreases as the 
sample size increases which indicates the asymptotic consistency property for the MLEs. It can be 
observed that as the reliability parameter R increases, the biases and MSEs for both R s x and R r x 

decreases. It is also observed that for small values of the scale parameter X, < A < 1 (Tables 1 and 2), 
estimation of R using both SRS and RSS approaches is better than estimation when A > 1 (Table3). 

Comparing the two approaches, the biases and MSEs of R r x is always better than those of R s ^ 

which can be noted from the efficiency of R r x with respect to R s x . Finally, the efficiency of R r x is 

greater than one and increases as the sample sizes increase. 

Table 1: MLE estimation of R when different values of the parameters a and /? 
and when the scale parameter is known ( A — 0.5 ) and r x — r = 5. 

R (n,m) m x ,m y Bias(R sX ) Bias(R rX ) MSE( R sA ) MSE(R rX ) eff(R rA ) 

0.250 (10,10) (2j2) 0.0699 0.0546 1.6481 1.4317 1.1511 

(10,15) (2,3) 0.0635 0.0496 1.2180 0.9648 1.2626 

(10,25) (2,5) 0.0585 0.0477 0.9625 0.7376 1.3056 

2358 | Vol. 6, Issue 6, pp. 2354-2362 



International Journal of Advances in Engineering & Technology, Jan. 2014. 

©IJAET ISSN: 22311963 





(15,25) 


(3,5) 


0.0551 


0.0441 


0.7336 


0.5507 


1.3319 




(25,25) 


(5,5) 


0.0512 


0.0398 


0.6122 


0.4154 


1.4738 


0.333 


(10,10) 


(2,2) 


0.0543 


0.0354 


1.5443 


1.3605 


1.1353 




(10,15) 


(2,3) 


0.0429 


0.0291 


1.11885 


0.9167 


1.2204 




(10,25) 


(2,5) 


0.0409 


0.0267 


0.8958 


0.7006 


1.2786 




(15,25) 


(3,5) 


0.0389 


0.0249 


0.6817 


0.5235 


1.3024 




(25,25) 


(5,5) 


0.0325 


0.0221 


0.5704 


0.3947 


1.4453 


0.500 


(10,10) 


(2,2) 


0.0423 


0.0305 


1.4764 


1.2079 


1.2223 




(10,15) 


(2,3) 


0.0358 


0.0258 


1.0912 


0.8137 


1.3405 




(10,25) 


(2,5) 


0.0299 


0.0225 


0.8623 


0.6219 


1.3864 




(15,25) 


(3,5) 


0.0271 


0.0203 


0.6572 


0.4646 


1.4143 




(25,25) 


(5,5) 


0.0243 


0.0196 


0.5484 


0.3504 


1.5650 


0.667 


(10,10) 


(2,2) 


0.0391 


0.0254 


1.3985 


1.2067 


1.1589 




(10,15) 


(2,3) 


0.0308 


0.0209 


1.0023 


0.7891 


1.2702 




(10,25) 


(2,5) 


0.0224 


0.0193 


0.7787 


0.5982 


1.3018 




(15,25) 


(3,5) 


0.0209 


0.0169 


0.5815 


0.4419 


1.3155 




(25,25) 


(5,5) 


0.0183 


0.0138 


0.4741 


0.3297 


1.4381 


0.850 


(10,10) 


(2,2) 


0.0269 


0.0175 


0.9367 


0.6756 


1.3865 




(10,15) 


(2,3) 


0.0213 


0.0145 


0.6786 


0.4552 


1.4908 




(10,25) 


(2,5) 


0.0203 


0.0126 


0.5433 


0.3479 


1.5617 




(15,25) 


(3,5) 


0.0173 


0.0108 


0.4135 


0.2598 


1.5916 




(25,25) 


(5,5) 


0.0154 


0.0083 


0.3461 


0.1962 


1.7640 


Table 2: MLE estimation of R when different values of the parameters a 


and ,3 




and when the scale parameter is known ( A — 1 ) and r x = r y = 5. 




R 


(n,m) 


m^m y 


Bias( R s X ) 


Bias( R rX ) 


MSE(R sX 


) MSE(R rX ) 


eff( R r,A ) 


0.250 


(10,10) 


(2,2) 


0.0572 


0.0404 


1.3494 


1.0606 


1.2722 




(10,15) 


(2,3) 


0.0510 


0.0367 


0.9775 


0.7147 


1.3676 




(10,25) 


(2,5) 


0.0429 


0.0320 


0.7059 


0.4944 


1.4278 




(15,25) 


(3,5) 


0.0412 


0.0299 


0.5489 


0.3729 


1.4722 




(25,25) 


(5,5) 


0.0361 


0.0230 


0.4314 


0.2397 


1.8001 


0.333 


(10,10) 


(2,2) 


0.0445 


0.0262 


1.2644 


1.0079 


1.2545 




(10,15) 


(2,3) 


0.0368 


0.0216 


0.9609 


0.6791 


1.4149 




(10,25) 


(2,5) 


0.0300 


0.0179 


0.6570 


0.4696 


1.3990 




(15,25) 


(3,5) 


0.0291 


0.0169 


0.5101 


0.3544 


1.4391 




(25,25) 


(5,5) 


0.0229 


0.0128 


0.4020 


0.2277 


1.7651 


0.500 


(10,10) 


(2,2) 


0.0346 


0.0226 


1.2088 


0.8948 


1.3508 




(10,15) 


(2,3) 


0.0287 


0.0191 


0.8757 


0.6028 


1.4527 




(10,25) 


(2,5) 


0.0219 


0.0151 


0.6325 


0.4169 


1.5171 




(15,25) 


(3,5) 


0.0203 


0.0137 


0.4918 


0.3146 


1.5633 




(25,25) 


(5,5) 


0.0171 


0.0113 


0.3865 


0.2022 


1.9116 


0.667 


(10,10) 


(2,2) 


0.0320 


0.0188 


1.1450 


0.8939 


1.2808 




(10,15) 


(2,3) 


0.0247 


0.0155 


0.8044 


0.5846 


1.3760 




(10,25) 


(2,5) 


0.0164 


0.0129 


0.5711 


0.4010 


1.4243 




(15,25) 


(3,5) 


0.0156 


0.0114 


0.4351 


0.2992 


1.4543 




(25,25) 


(5,5) 


0.0129 


0.0080 


0.3341 


0.1902 


1.7563 


0.850 


(10,10) 


(2,2) 


0.0220 


0.0130 


0.7669 


0.5005 


1.5323 




(10,15) 


(2,3) 


0.0171 


0.0107 


0.5446 


0.3372 


1.6149 




(10,25) 


(2,5) 


0.0149 


0.0084 


0.3985 


0.2332 


1.7087 




(15,25) 


(3,5) 


0.0129 


0.0073 


0.3094 


0.1759 


1.7590 




(25,25) 


(5,5) 


0.0109 


0.0048 


0.2439 


0.1132 


2.1546 


Table 3: MLE estimation of R when different values of the parameters a and /? 






and when the scale parameter is known ( X = 


= 3) and r x = r y 


= 5. 


R 


(w,m) 


m x ,m y 


Bias(R sX ) 


Bias( R rX ) 


MSE(R sX 


) MSE(R rX ) 


eff(RrA) 


0.250 


(10,10) 


(2,2) 


0.4673 


0.2989 


3.0030 


2.3605 


1.2722 




(10,15) 


(2,3) 


0.3842 


0.2457 


2.0081 


1.4393 


1.3952 




(10,25) 


(2,5) 


0.3202 


0.2138 


1.4359 


0.9957 


1.4421 




(15,25) 


(3,5) 


0.2729 


0.1788 


0.9903 


0.6726 


1.4722 




(25,25) 


(5,5) 


0.2295 


0.1460 


0.7477 


0.4591 


1.6288 


0.333 


(10,10) 


(2,2) 


0.3630 


0.1938 


2.8139 


2.2431 


1.2545 




(10,15) 


(2,3) 


0.2595 


0.1441 


1.8447 


1.3676 


1.3489 




(10,25) 


(2,5) 


0.2239 


0.1197 


1.3364 


0.9457 


1.4131 
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(15,25) 


(3,5) 


0.1927 


0.1010 


0.9202 


0.6394 


1.4391 


(25,25) 


(5,5) 


0.1457 


0.0811 


0.6967 


0.4362 


1.5971 


(10,10) 


(2,2) 


0.2828 


0.1670 


2.6902 


1.9915 


1.3508 


(10,15) 


(2,3) 


0.2166 


0.1278 


1.7991 


1.2139 


1.4821 


(10,25) 


(2,5) 


0.1637 


0.1008 


1.2864 


0.8395 


1.5324 


(15,25) 


(3,5) 


0.1342 


0.0823 


0.8871 


0.5675 


1.5633 


(25,25) 


(5,5) 


0.1089 


0.0719 


0.6698 


0.3873 


1.7297 


(10,10) 


(2,2) 


0.2614 


0.1390 


2.5482 


1.9895 


1.2808 


(10,15) 


(2,3) 


0.1863 


0.1035 


1.6525 


1.1772 


1.4038 


(10,25) 


(2,5) 


0.1226 


0.0865 


1.1617 


0.8075 


1.4386 


(15,25) 


(3,5) 


0.1035 


0.0685 


0.7849 


0.5397 


1.4543 


(25,25) 


(5,5) 


0.0820 


0.0506 


0.5791 


0.3644 


1.5892 


(10,10) 


(2,2) 


0.1799 


0.0958 


1.7068 


1.1139 


1.5323 


(10,15) 


(2,3) 


0.1289 


0.0718 


1.1188 


0.6791 


1.6476 


(10,25) 


(2,5) 


0.1111 


0.0565 


0.8105 


0.4696 


1.7259 


(15,25) 


(3,5) 


0.0857 


0.0438 


0.5582 


0.3173 


1.7590 


(25,25) 


(5,5) 


0.0690 


0.0305 


0.4227 


0.2168 


1.9495 



VI. Conclusions 

In this paper, we have addressed the problem of estimating R = P(Y < X) for the generalized inverted 
exponential distribution using both simple random sampling SRS and ranked set sampling RSS 
approaches. Comparing the performance of estimators of R, it is observed that: 

1 . The biases of all estimators is very small. 

2. The estimators based on RSS dominate those based on SRS. 

3. Relative to biases and MSEs, the ML estimator of R based on RSS is better than that based on 
SRS for different a, P and X parameter values 

4. Similarly, the Bayes estimator of R based on RSS is better than that based on SRS for 
different a, P and X parameter values 

5. Bayes estimator of R based on both approaches dominate the corresponding ML estimators. 

6. The MSEs of estimators made using RSS is always less than those made by SRS. 

7. Though out the simulation studies the estimation of R using the RSS approach is more 
efficient than that made using the SRS approach. 

8. The efficiency of all estimators of based on RSS is decreasing as the set sizes is increasing. 

9. This study revealed that the estimator of P(Y < X) based on RSS technique is more efficient 
than the corresponding estimator based on SRS technique. 



VII. Future Work 

Many devices cannot function at high as well as low temperatures and a person's blood pressure has 
two limits (systolic and diastolic pressures). This case can be represented as the reliability parameter 
7?3=P(Xi<Y<X2), where Y is the stress that is limited by two strengths Xi , and X2- In the near future 
estimation R3 is considered for several lifetime distribution with real data application and using 
several sampling approaches. 
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